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differential  algebraic  equations.  His  thesis  is  titled  Constraint  Preserving  Integrators  for 
Differential  Algebraic  Equations  and  will  ccxnplete  in  the  fall  of  1993.  Mr.  Moore  will 
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9.  REPORT  OF  INVENTIONS  (BY  TITLE  ONLY): 

There  were  no  inventions  daring  this  period. 

BRIEF  OUTUNE  OF  RESEAROB  FINDINGS 

Differential  algebraic  equations  (DAEs)  are  systems  cl  differential  and  algebraic  equar 
tions  which  have  not  been  reworked  into  an  explicit  format.  Many  engineering  problems  are 
initially  modeled  as  DAEs.  Reduction  to  explicit  models  can  be  parameter  dependent,  time 
consuming,  and  require  model  simplification.  Small  design  changes  can  require  a  conq>lete 
rederivatkm  of  an  explicit  model.  This  project  is  to  develop  numerical  methods  for  working 
directly  with  the  original  implicit  eqnatkms. 

Previously  developed  numerical  methods  for  DAEs  such  as  backward  differentiation 
(BDF)  and  implicit  Runge-Kutta  (IRK)  methods  experience  numerical  difficulties  with 
the  higher  index  problems  that  arise  in  such  areas  as  constrained  mechanics,  for  example 
with  models  for  vehicular  systems.  Even  when  these  numerical  methods  can  be  used  it  is 
important,  and  often  difficult,  to  determine  sufficiently  accurate  initial  conditions  for  all 
variables. 

In  1988,  one  the  principal  investigators  outlined  a  more  general  numerical  method  for 
nonlinear  DAEs.  While  computationally  mote  intensive  than  BDF  and  IRK  methods,  it 
had  several  potential  advantages.  It  will  work  when  the  other  methods  do  not.  It  can  also 
be  used  to  compute  consistent  initial  conditions  fw  use  with  either  this  method  or  for  BDF 
and  IRK  integratcns  when  they  are  apprc^riate.  This  general  method  invcdved  symbolic 
differentiations  of  the  original  equations  followed  by  a  singular  Newt(»’s  method  which  was 
carried  out  numerically. 

One  emphasis  of  this  project  has  been  to  develop  the  theoretical  basis  of  the  general 
method.  Tliis  has  proven  to  be  highly  n(»trivial  but  substantial  progress  has  been  made. 
In  [3]  we  showed  that  nonlinear  index  three  Hessenberg  systems  can  solved  by  this  method. 
This  result  is  important  since  many  mechanics  and  cmtrol  problems  have  this  form. 

A  closely  related  aspect  of  our  research  has  been  to  carefully  examine  the  information 
present  in  the  derivative  arrays  which  occur  in  the  general  method  [1,3,8,9,11,15,18.19]. 
This  information  will  eventually  be  used  to  develop  nonlinear  boundary  value  (BVP)  codes 
along  the  lines  suggested  by  CHark  and  Petzold.  However,  we  have  initially  applied  our 
results  to  determining  ccmtrol  thewetic  information  f<»  implicit  contrrd  problems  [1].  As 
shown  in  our  earlier  work,  this  will  have  direct  explication  to  such  problems  as  tracking  or 
path  control  for  mechanical  systems.  These  results  are  also  relevant  for  classical  methods 
since  porti<ms  of  the  derivative  arre^  are  often  used  to  reformulate  higher  index  DAEs. 

In  other  related  research,  we  have  devel<x^  an  algorithm  for  computing  the  state  space 
decomposition  of  linear  time  varying  Hessenberg  systems  of  arbitrary  index  [5,10].  This 
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algorithm  does  not  involve  derivative  array  computationa  but  does  require  8ynd>olic  differ- 
entiation  of  the  coefficients  and  input  functions.  In  principle  it  can  be  used  to  c<nnpute 
exactly  ccwsistent  initial  conditions  for  linear  Heasenberg  systems,  and  ccmsequently  to  de¬ 
termine  c<»nplementary  omditkms  for  the  unique  scdvability  of  two-point  boundary  value 
problems  for  these  systems.  An  analogous  algorithm  can  be  used  to  analsrse  the  errors  in 
BDF  methods  f<»  arbitrary  index  Heasenberg  DAEs  to  establish  convergence  under  suitable 
smoothness  aasuiiq>tions  on  the  starting  values.  These  are  the  first  results  are  linear  Hes- 
senberg  systems  without  any  restrictions  on  the  index  and  resolve  several  open  questi<»s  in 
the  literature. 

A  sec<»d  thrust  of  our  research  has  been  to  develop  numerical  methods  and  guidelines 
for  their  use.  Mr.  Mocne,  the  research  assistant,  is  currently  developing  a  software  package 
that  wiU  carry  out  the  entire  general  procedure  for  both  integrating  DAEs  and  for  finding 
initial  conditions.  The  user  will  need  only  to  write  down  the  equatkms  and  specified  desired 
quantities.*  The  software  will  then  symbolically  differentiate  the  equations  several  times, 
generate  Jacobians  of  the  eorreq>onding  derivative  arrays,  generate  FORTRAN  code  fi» 
these  arrays  and  Jacobians,  and  finally  carry  out  the  numerical  calculations  in  FORTRAN. 
One  advantage  of  this  software  is  that  the  symlxdic  part  wiU  need  to  be  done  only  once. 
Design  and  model  parameters  can  then  be  varied  in  the  numerical  routines.  [6]  is  the  first 
report  on  this  prototype  code.  Numerical  results  reported  in  [6]  on  the  shuttle  reentry 
problem  of  Brenan  show  that  contrary  to  some  prior  expectations,  the  growth  of  symbolic 
expressions  under  repeated  differentiations  will  not  be  such  a  maj<»  factor  with  the  method 
thanks  to  the  new  optimise  feature  in  the  fortran  command  in  MAPLE  V.  Rather,  the 
dominant  factor  will  usually  be  the  dimension  of  the  problem.  Indications  are  that  the 
prototype  code  will  be  able  to  handle  moderate  sised  nonlinear  DAEs  (10-20  variables)  of 
higher  index  (1-6)  in  reasonable  time. 

In  many  mechanical  and  control  systems  their  ate  delays  present.  There  has  been 
essentially  no  prior  wcnk  none  on  either  tlM  numerical  or  analytic  s<dution  of  differential- 
delay  DAEs.  Initial  results  on  this  problem  are  in  the  invited  [4,13]. 

Regardless  of  the  integrator  being  used  for  general  DAEs  some  type  of  differentiatimi 
is  needed.  While  we  have  shown  that  this  can  be  d(»e  f<Mr  many  problems  with  reascmable 
speed,  H  is  hi^^  desirable  to  be  able  to  speedup  this  part  of  the  computation.  This 
will  permit  faster  integration  and  the  consideration  of  significantly  larger  problems.  One 
software  techndogy  that  looks  promising  in  this  regard  is  automatic  differentiation.  In 
coqjunction  with  Mr.  Zhmig,  a  Ph.D.  student  of  one  of  the  Pi’s,  we  investigated  the  use  of 
autcnnatic  differentiation  with  DAEs.  [12]  repmts  on  our  experience  with  the  ADOL-C  code 
of  Dr.  Griewank  of  Argonne  National  Laboratory.  We  found  that  automatic  differentiatimi 
can  result  in  a  significant  speed  up  of  the  ccmiputation,  and  a  reduction  in  storage,  under 
some  operating  ccmditions.  In  other  scenarios  automatic  differentiatimi  was  actually  more 
time  cmisuming. 
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Tlie  general  approach  described  in  previous  reports  has  several  advantages.  However, 
it  has  the  disadvantage  that  it  does  not  rein^Kiee  the  explicit  and  implicit  constraints  that 
are  present  in  a  DAE.  Thus  there  can  be  drift  off  the  cmistraint  manifcdd.  This  is  c^n 
undesirable,  particularly  in  siinulatkma  of  vehicular  and  other  systems.  Mr.  Moore  has 
examined  the  use  of  the  ideas  from  the  general  ^proach  to  design  an  integrator  which  will 
work  for  unstructured  higher  index  DAEs  and  also  omserve  constraints.  In  [18]  we  outline 
the  first  general  higher  index  DAE  integrat<»  which  preserves  constraints.  The  method  is 
perfuming  well  in  tests  and  has  already  integrated  several  DAEs  that  other  methods  cannot 
integrate. 

Most  of  the  prior  work  on  DAEs  has  assumed  that  the  system  was  stdvable,  that  is 
that  there  is  a  well  defined  family  ci  solutkms.  However,  in  applicatkms  there  can  occur 
situations  where  the  system  is  not  solvable.  In  mechanics  this  occurs,  fu  example,  when 
the  object  is  in  a  singular  configuratkm  and  there  is  a  sudden  change  in  the  degrees  of 
freedom.  Fu  a  cunplex  system  it  may  not  be  obvious  when  this  occurs.  Recently  there 
has  been  sune  work  in  the  mathematical  literature  by  Rheinboldt,  Rabier,  and  Reich  cm 
proving  solvability  for  nonlinear  DAEs.  This  wuk,  while  important,  is  based  on  a  recursive 
nonlinear  verskm  cff  the  arguments  used  for  the  linear  time  invariant  case.  It  can  be  very 
difficult  to  actually  carry  out  these  procedures.  In  [8]  we  go  a  kmg  way  toward  rectifying 
this  problem.  We  give  direct  characterisations  of  solvability  that  are  not  recursive  in  nature. 
These  characterisations  are  under  weaker  assumptions  than  any  other  characterisations  in 
the  literature.  We  also  show  how  cme  can  actually  prove  solvability  on  a  computer  using 
a  mixture  oi  readily  available  symbdUc  and  numeric  software.  An  example  is  given  frcun 
mechanics  ccmsisting  of  the  path  control  of  a  robotic  arm  with  a  flexible  joint.  The  standard 
graph  theoretic  ^proaches  give  an  incorrect  answer  on  this  problem  but  our  approach  easfly 
handles  it. 

Many  problems,  particularly  in  mechanics,  possess  a  structure  which  is  known  as  the 
Hessenberg  structure.  Index  two  and  index  three  versions  of  this  structure  have  been  utilised 
throu^out  the  numerical  literature  on  DAEs.  However,  mechanical  and  other  systems 
can  have  even  higher  index.  In  [5]  we  have  established,  for  the  first  time,  the  necessary 
mathematical  machinery  f<v  analysing  problems  with  Hessenberg  systems  of  arbitrary  index. 
In  [10]  we  iH;>ply  this  mathematical  machinery  to  analysing  the  behavior  of  the  numerical 
BDF  (Backward  Differentiation  Formulas)  methods  on  these  systems.  These  results  provide 
the  definitive  result  on  these  systems  which  have  been  studied  for  the  last  decade  in  special 
8ituati<»8.  The  results  should  also  prove  useful  for  other  problems  such  as  the  c(Hitrol  of 
Hessenberg  type  systems. 

Dr.  Marek  Rakowski  visited  NCSU.  Jointly  he  and  one  of  the  Pis  examined  whether  his 
results  on  rational  matrix  functions  could  be  useful  for  DAEs.  In  [7]  they  showed  that  this 
theory  could  be  used  to  find  tsplicU  completims  of  linear  time  varying  DAEs  with  rational 
matrix  coefficients.  While  the  general  usefulness  of  this  work  is  not  clear,  it  will  have  at 
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least  omt  inqixMrtaiit  benefit.  It  will  naake  poauble  the  derivati<»  of  exact  soluticHu  for  very 
c(mq>lex  Imear  time  varying  DAEs.  This  will  be  very  useful  for  deriving  test  problems  for 
testing  more  general  numerical  approaches  <». 

Building  on  the  previous  papers  two  additional  important  matbematical  contributkms 
have  been  made  to  the  understanding  of  DAEs  in  the  last  year.  There  has  been  a  perception 
in  the  Hterature  that  most  at  the  concepts  o(  index  were  equivalent.  Spurred  by  the  need 
to  consider  more  cmnplex  systems  we  investigated  this  carefully  in  [15].  There  we  showed 
that  the  different  types  of  index  could  be  very  different  in  value  for  fully  implicit  DAEs. 
We  also  provided  an  unifying  thecwy  and  related  these  ideas  to  the  numerical  results  in  [8]. 
In  addition,  it  is  standard  practice  in  engineering  to  linearise  equations.  Little  work  had 
been  done  rm  the  linearisations  of  hi^er  index  DAEs.  In  [19]  we  showed  that  the  usual 
time  invariant  linearisations  about  trajectories  were  inapprc^riate  but  that  time  varying 
linearisations  could  be  ^prr^riate  for  srune  applications. 

In  order  tot  the  research  <»  DAEs  to  he  useful  it  is  imp<»tant  ftx  the  results  to  be  broui^t 
to  the  attention  of  scientists  and  engineers,  many  of  wh(»n  do  not  read  the  mathematical 
journals  where  the  results  are  inhialfy  published.  ^,14,17],  which  are  primarily  e^qxisitory, 
are  part  of  our  effort  to  inform  this  larger  audience  of  the  current  successes,  and  long  term 
prospects,  oi  the  current  research  on  implicit  problems  and  their  numerical  stdution.  There 
is  a  misccmception  in  some  scientific  cmnmunities  that  higher  index  DAEs  are  not  common. 
In  [20]  we  show  that  they  can  frequently  occur,  particularly  when  more  c<Mnplex  models 
incorporating  actnaUw  and  other  cmitrol  affects  are  included. 

POTENTIAL  TECHNOLOGICAL  APPLICATIONS 

The  early  work  on  index  one  DA!'  .  lead  to  the  development  of  such  numerical  production 
codes  as  LSODI  and  DASSL,  and  more  recent  codes  such  as  DASSLRT  and  RADUAIIa. 
These  codes  are  being  used  daily  at  lab<»at<mes  for  problems  from  chemical  engineering  to 
open  lo<^  mechanical  system  simulatkm.  It  is  expected  that  the  research  of  this  project, 
combined  with  that  going  <»  elsewhere  in  the  United  States  and  in  Europe,  will  lead  to  the 
development  of  production  nonlinear  hi^er  index  codes  within  the  next  few  years.  These 
codes  will  have  immediate  application  in  several  areas.  One  such  area  is  the  design  and 
simulati<»  of  cmnplex  mnhihody  dynamical  systems  such  as  moving  multihody  vehicular 
systems  and  telerohots. 

CONTACTS  WITH  ARMY  LABORATORY  PERSONNEL 

The  Pis  made  a  two  day  visit  in  1991  to  the  US  Army  Tsnk-Autmiiotive  Command 
(TACOM)  May  1  &  2.  During  thu  visit  they  talked  with  Dr.  R.  A.  Wehage,  Dr.  N. 
A.  Rose,  and  Mr.  J.  L.  Overman.  The  visit  was  to  familiarise  the  Pis  with  the  pers<»nel 
and  research  interests  of  this  group  at  TACOM.  Particular  topics  of  discussion  were  the 
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munerical  nmnlation  of  mechanical  systems  (Wehai^e),  nonlinear  control  of  the  Stewart 
Platfonn  Sinnilatm  (Overhoh),  and  a  numerical  problem  arising  in  optics  (Rose). 

Dt.  Rose  had  a  problem  that  he  thought  might  be  handled  better  symbolically.  One 
oi  the  Pi’s  (S.  CampbeD)  subsequently  held  two  telephone  conversations  with  Dr.  Rose 
and  ran  the  synbolk  calculatkms  using  the  symlxdic  language  MAPLE  oa  a  workstation 
at  NCSU.  The  expressions  turned  out  to  be  so  complicated,  however,  that  it  was  decided 
that  it  would  be  better  to  just  evaluate  them  by  standard  iq>proaches. 

Since  the  numerical  solution  of  DAEs  is  of  direct  interest  to  the  Axvay  for  such  pur¬ 
poses  as  vehicular  simulatkm,  a  timely  transfer  of  results  and  problems  between  labwat<»y 
personnel  and  university  researchers  is  implant.  Drs.  Clark  and  Caiiq>bell,  in  coordinar 
tion  with  Dt.  Linda  Petsold  ci  the  Univetmty  of  Minnesota,  <»ganized  an  ARC  workshop 
on  Modelling  and  Analysis  for  Mechanical  Systems  and  Algorithms  for  Real  Time  Simu¬ 
lation.  The  wcwkshcq)  was  held  November  5-7,  1992,  at  the  Brownstone  Hotel  in  Raleigh, 
NC.  Thete-were  several  Army  Lahoratcwy  personnel  involved  in  this  workshop  as  well  as 
University  persMinel  and  sune  researchers  from  laboratories  in  Europe.  All  indications  are 
that  the  w<»ksh<^  was  viewed  as  highly  successful. 

In  addition  Dr.  Clark  will  be  editing  a  special  issue  of  the  Journal  of  Mechanics  of 
Structures  and  Machines  which  will  consist  of  piq>ers  presented  at,  or  arising  out  of,  this 
workshop.  Dr.  Roger  Wehage  of  TACOM  may  be  assisting  Dr.  Clark  in  the  editing. 

ABSTRACTS  OF  RESEARCH  MANUSCRIPTS 

1.  S.  L.  Caiiq>bell,  N.  K.  Nichcds,  and  W.  J.  Terrell,  DuaUtf,  ohtemhilHji,  and  eonirol- 
UitUijf  for  Hnear  time  varfing  iescrifior  systems 

A  characterisatkm  of  observability  for  linear  time  varying  descriptor  systems  E{t)x'{t) 
-|-F(t)c(<)  s  B(<)u(<),  y(t)  =  C(t)x(<),  was  recently  developed.  Neither  E  nor  C  were 
required  to  have  cmistant  rank.  This  paper  defines  a  dual  system,  and  a  type  of  c<»i- 
trollability  so  that  observability  of  the  original  system  is  equivalent  to  cmitrcdlability 
of  the  dual  system.  Criteria  for  observability  and  controllability  are  given  in  terms  of 
arrays  of  derivatives  of  the  <»iginal  coefficients.  In  addition,  the  duality  results  cf  this 
p^>er  lead  to  an  improvement  on  a  previotu  frmdamental  structure  result  for  scdvable 
systems  of  the  ftwm  £^(t)x'(t)  +  F'(t)x(<)  =  /(<). 

2.  S.  L.  Campbell,  Descriptor  Spstems  ta  Ike  90 ’s 

After  briefly  summarising  the  ^owth  in  the  theory  of  deaerator,  or  differential  alge¬ 
braic  equatkms  (DAEs),  over  the  last  two  decades,  the  current  challenges  and  potential 
successes  will  be  discussed.  This  paper  was  the  introductory  paper  few  a  special  seaskm 
<»  DAEs  at  the  29th  IEEE  Conference  on  Decision  and  Cmitrol,  Dec.  1990. 
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3.  S.  L.  Campbell,  Ltsst  cotnpUUom  of  nonlinear  index  three  Hesaenherg  DAE$ 

Many  control  problems  are  most  easily  initially  modeled  as  a  nonlinear  implicit  system 
of  (Merential  and  algebraic  equations  (DAEs),  =  0.  Recently  a 

numerical  metbod  has  been  prc^osed  that  can  be  used  when  classical  methods  cannot 
be  used.  This  invited  p^r  e<»tinues  the  examination  of  bow  general  that  metbod  k. 
In  particular,  it  k  shown  that  nonlinear  index  three  Heasenberg  DAEs  can  be  solved 
by  tbk  method. 

4.  S.  L.  Campbell,  2-D  (Differential  Delap)  Implidi  apatemty 

Tbk  invited  p^>er  discusses  what  k  known  about  linear  time  invariant  differential- 
delay  descriptor  systems.  The  literature  k  surveyed  and  several  observatkms  are 
made. 

5.  K.  D.  Clark,  Decompoaition  of  Heaaenherg  DAE  apatema  to  state  apace  form 

An  algorithm  k  given  for  syiiy[>o]ically  decoupling  the  solutions  of  a  linear,  time  depen¬ 
dent  DAE,  Ez*  s  A(i)z  +  f(t)y  z{t)  €  R'  in  Hessenberg  form  into  state  and  algebraic 
coiiq>onent8.  The  state  variables  are  the  solutions  to  an  ordinary  differential  equation 
with  initial  conditions  restricted  to  a  subspace  ai  R*  while  the  algebraic  components 
are  linear  ffinctions  of  the  state  variables  which  involve  derivatives  of  the  coefficients 
and  input  fonctimis  up  to  order  r- 1  where  r  k  the  index  of  the  system.  The  iiiq>licar 
tions  of  the  algOTithm  for  computing  ccmsktent  initial  conditions,  for  certain  singular 
q>timal  ccmtrcd  problems,  and  for  nun^rical  solutions  are  briefly  dkcussed. 

6.  S.  L.  Campbell  and  E.  Moore,  Progreaa  on  a  general  numerical  method  for  nonlinear 
higher  index  DAEa 

A  method  has  been  prq[>06ed  for  numerically  solving  lower  dimensicmal,  nonlinear, 
higher  index  differential  algebraic  equations  for  which  more  classical  methods  such 
as  backward  differentiation  or  iiiq>licit  Runge-Kutta  may  not  be  appropriate.  Thk 
method  k  baaed  on  solving  nonlinear  DAE  derivative  arrays.  Thk  paper  discusses 
prc^ess  <m  the  implementati<»  of  thk  method,  resolves  some  of  the  issues  involved, 
and  Ikts  s<xne  remaining  problems.  Computational  experience  showing  that  the  «q>- 
proach  should  prove  practical  for  many  ^phcations  k  presented. 

7.  S.  L.  Campbell  and  M.  Rakowski,  Explicit  formulae  for  compkiiona  of  linear  time 
varping  aingular  apatema  of  differential  eguationa 

Completions  of  linear  time  varying  singular  systems  of  the  form  = 

f{t)  are  explicitly  cmnputed  using  recent  results  on  rational  matrix  fimctions.  The 
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algwithm  and  the  theory  behind  it  are  carefully  described.  Computational  issues  are 
discussed. 

8.  S.  L.  Campbell  and  E.  Griepentrog  SchthiiUjf  of  gatent  iiffereniiol  algehnie  egua- 
iious 

In  tlM  last  few  years  there  has  been  ccmsiderable  research  mi  differential  algebraic 
equations  (DAEs)  /(<,c,x')  s  0  where  fx^  is  identically  singular.  Most  of  this  effmi 
has  focused  on  computing  a  solutimi  that  is  assumed  to  exist.  Mme  recently  there 
have  been  existence  results  developed  using  differential  geometry  in  a  recursive  manner 
that  fi^ows  a  pattern  similar  to  that  used  in  the  linear  time  invariant  case.  In  this 
paper  we  shall  establish  the  existence  of  scdutimis  using  more  direct  characterixatkms. 
Verificatimi  of  these  cmiditions  using  readily  available  numerical  and  symbolic  software 
is  discussed.  An  example  from  robotics  where  classical  graph  theoretical  approaches 
give  an  incewrect  answer  is  worked  to  illustrate  the  usefulness  df  the  characterization. 

9.  S.  L.  Campbell,  Unigutness  of  CompUiions  for  linear  time  varying  iiffereniial  alge- 
hraie  equations 

The  extent  to  which  a  cmnpletion  x*  s=  Git)x + Ri(t)f^'\t)  of  a  linear  time  vary¬ 
ing  differential  equaticm  E(t)x'{t)  +  F(t)x(t)  s  /(<)  is  unique  is  carefully  examined. 
The  implications  for  numerical  methods  for  solving  DAEs  based  on  differentiated 
equatkms  are  discussed. 

10.  K.  D.  Clark,  A  framework  for  arbitrary  index  differential-  algebraic  egnations  in  Hes- 
senberg  form 

This  piqier  presents  an  algebraic  framewmk,  related  to  a  state  space  decomposition 
algorithm,  which  can  be  used  to  analyze  the  cmivergence  behavior  of  backward  dif¬ 
ferentiation  formulas  for  the  solution  of  differential-algebraic  equatims  in  Hessenberg 
fmm.  The  system  cmisidered  is  assumed  to  be  linear  time  varying  but  is  allowed  to 
have  arbitrary  fixed  index.  In  this  prehminary  repmt,  the  index  two  and  three  cases 
are  discussed  in  detail  to  illustrate  both  the  usefulness  of  the  iqiproach  deveh^ed, 
and  also  the  increase  in  complexity  of  the  analysb  as  the  index  increases.  However, 
the  growth  of  cmnplexity  is  regular  in  such  as  manner  as  to  allow  for  the  general 
case  analysis.  In  the  final  versimi  of  this  paper,  the  analysis  will  be  extended  to  the 
arbitrary  index  linear  case.  The  nonlinear  case  will  be  given  in  a  subsequent  paper. 

11.  S.  L.  CampbeU,  Least  squares  completions  for  nonlinear  differential  algebraic  equations 
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A  method  has  been  prt^oeed  for  numerically  solving  lower  dimensicmal,  nonlinear, 
hitter  index  differential  algebraic  equations  f<x  which  mme  classical  methods  such 
as  backward  differentiation  or  implicit  Runge-Kutta  may  not  be  appropriate.  This 
method  is  baaed  on  solving  n<mlinear  DAE  derivative  arrays  using  nonlinear  singular 
least  squares  methods.  The  theoretical  foundations,  generality,  and  limitations  of  this 
iq>proach  remain  to  be  determined.  This  p^er  carefully  examines  several  key  aspects 
ci  this  approach.  The  emphasis  is  on  general  results  rather  than  specific  results  based 
<»  the  structure  of  various  ^plications. 

12.  S.  L.  Campbell,  E.  Moore,  and  Y.  Zhcmg,  Uiilizaiion  of  avtomatic  Hfftreniiaiion  tn 
control  nlforiikma 

Symbi^c  languages  are  increasingly  being  used  in  the  analysis  and  implementation  of 
control  algorithms.  Many  of  these  control  procedures  invcdve  scmie  typo  ct  differentia¬ 
tion  at  Jacobian  formulation.  Automatic  differentiation  provides  an  alternative  means 
of  c<xnputing  this  infcwmation  which  is  rarely  considered  in  the  contrd  literature.  This 
piq>er  will  discuss  the  use  of  automatic  differentiation  in  c(»trol  problems.  Numerical 
results  comparing  automatic  and  symbolic  ^proaches  on  a  nonlinear  control  system 
will  be  given. 

13.  S.  L.  Campbell,  Nonrcgnlar  dcseripior  sjfatcms  with  dekya 

This  paper  examines  descriptor  systems  with  delays  which  are  in  the  form  Bx'(t)  + 
Dx{t) + Cx(t  —  or)  s  Eu{t).  Previous  papers  have  examined  the  case  when  XB  +  Dio 
a  regular  pencU.  In  this  paper  we  examine  the  case  when  XB  +  Dan  singular  pencil 
but  det(AB  +  D  +  wC)  ^  0.  Behavior  not  previously  reported  for  descriptor  systems 
is  observed.  Several  basic  issues  are  discussed. 

14.  S.  L.  Campbell,  A  survey  of  time  varying  and  nonlinear  descriptor  control  systems 

The  mathematical  results  on  the  structure  of  matrix  pencils  are  well  known  in  the 
control  community.  This  theory  has  been  heavily  exploited  by  control  theorists  work¬ 
ing  with  linear  time  invariant  control  problems  in  descriptew  form.  The  structure  and 
pr<^erties  of  time  varying  and  nonlinear  descriptor  systems  are  less  well  known  in 
the  systems  and  control  community  since  most  of  these  results  has  been  developed 
in  recent  years  and  published  outside  of  the  usual  control  literature  by  analysts  and 
numerical  analysts.  This  nonlinear  thewy  has  taken  on  a  very  different  appearance 
from  that  of  the  linear  time  invariant  theory.  This  survey  p^r  discusses  what  is 
known  about  nonlinear  and  time  varying  descriptor  systems  and  the  ways  in  which 
they  are  different  from  linear  time  invariant  descriptor  systems.  Both  analytical  re¬ 
sults  and  numerical  methods  will  be  mentioned.  Hopefully  this  will  give  the  listener 
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some  ideft  ci  what  took  are  available  for  attaddng  control  problems  and  what  types 
of  assumptions  may  be  necessary  in  (xder  to  apply  these  results. 

15.  S.  L.  Campbell  and  C.  W.  Gear,  The  index  of  general  nonlinear  DAEs 

In  the  last  few  years  there  has  been  considerable  research  on  differential  algebraic 
equations  (DAEs)  =  0  where  Fyt  is  identically  singular.  Most  of  this  effort 

has  focused  on  computing  a  solution  that  is  assumed  to  exist.  More  recently  there 
has  been  some  dkcusskm  of  solvability  of  DAEs.  There  has  historically  been  some 
inq>reckion  in  the  use  of  the  two  key  concepts  of  sdvability  and  index  for  DAEs.  The 
ccmsideration  of  increasingly  complex  fully  implicit  nonlinear  DAEk  makes  a  clear  and 
correct  development  necessary.  This  paper  will  try  to  clarify  several  points  concerning 
sdval.!Mty  and  the  index.  After  establishing  some  new  and  more  precise  terminology 
that  we  need,  some  inaccuracies  in  the  literature  wiU  be  corrected.  The  two  types 
of  indices  most  frequently  used,  the  differentiation  index  and  the  i>erturbation  index, 
are  defined  with  respect  to  scdnti<ms  of  unperturbed  problems.  Examples  are  given 
to  show  that  these  indices  can  be  very  different  for  the  same  problem.  We  define  new 
"maximum  indices,”  which  are  the  maxima  of  earlier  indices  in  a  neighborhood  of  the 
solution  over  a  set  of  perturbaticms  and  show  that  these  indices  are  simply  related  to 
each  other.  These  indices  are  also  related  to  an  index  defined  in  terms  of  Jacobians 
of  the  problem. 

16.  S.  L.  Campbell  and  Y.  Shtessel,  Sliding  mode  control  and  differeniitd  algebraic  egna- 
iiona 

Sliding  mode  controlkrs  can  be  advantageous  in  stabilization  because  of  increased 
robustness  under  small  disturbances.  For  some  systems,  implementation  of  a  sliding 
mode  controller  can  involve  the  solution  of  higher  index  differential  algebraic  equations 
(DAEs).  The  relationship  between  DAEs  and  sliding  mode  control  will  be  examined 
in  this  p^er.  A  voltage  stabilizer  example  is  carefully  considered. 

17.  S.  L.  Campbell,  Numerical  methods  for  nnstmetnred  higher  index  DAEs 

Many  physical  problems  are  initially  modeUed  as  systems  of  differential  and  algebraic 
equations  (DAEs).  The  widespread  occurrence  of  DAEs  has  lead  to  an  intensive  ex¬ 
amination  of  the  numerical  solution  of  DAE  initial  value  problems.  Several  numerical 
methods  have  been  developed  for  DAEs.  Sc«ne  of  these  methods  are  suited  either  for 
what  are  called  low  index  problems  or  systems  with  special  structure.  There  is  a  need 
for  numerical  methods  for  nonlinear  higher  index  DAEs  which  do  not  have  any  of  the 
obvious  structure  to  exploit.  This  paper  will  survey  what  is  currently  known  about 
general  DAEk  and  methods  for  their  numerical  solution. 
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18.  S.  L.  Caiiq>bell  and  E.  Moore,  Couatnini  prtaerving  inUgratora  for  general  nonlinear 
higher  index  DAEa 

In  the  last  few  years  there  has  been  e(»siderahle  research  cm  numerical  methods 
f<x  differential  algebraic  equaticms  (DAEs)  f{x',x,t)  =  0  where  is  identically 
singular.  The  index  provides  erne  measure  of  the  singularity  of  a  DAE.  Most  of  the 
numerical  analysis  literature  cm  DAEs  to  date  has  dealt  with  DAEs  with  indices  no 
larger  than  three.  Even  in  this  case,  the  systems  were  often  assumed  to  have  a  special 
structure.  Recently  a  numerical  method  was  pr(q>osed  that  could,  in  principle,  be 
used  to  integrate  general  unstructured  hi|^er  index  solvable  DAEs.  However,  that 
method  did  not  preserve  cemstramts.  This  p^er  will  discuss  a  modiffcation  oi  that 
approach  which  can  be  used  to  design  constraint  preserving  integrators  for  general 
nonlinear  higher  index  DAEs. 

19.  S.  L. 'Campbell,  Lineariaaiion  of  DAEa  along  irajedoriea 

Over  the  last  decade  there  has  been  a  considerable  amount  of  research  cm  numerical 
and  analytic  aspects  of  linear  and  nemlinear  differential  algebraic  equations  (DAEs) 
F(x*,x,t,u)  s  0.  Many  of  these  papers  have  either  cemsidered  linear  equations  or 
based  their  analysis  cm  linear  equatiems.  However,  until  very  recently  there  has  been 
little  rigorous  investigatiem  of  the  relationship  between  the  linearization  of  a  DAE 
and  the  criginal  equatiems.  In  this  paptt  we  carefully  examine  several  aspects  of  this 
relationship.  Positive  results  for  time  varying  linearization  and  counter  examples  for 
time  invariant  linearization  are  given. 

20.  S.  L.  Campbell,  High  index  differential  algebraic  egnaiiona 

In  the  last  few  years  there  has  been  ccmsiderable  research  on  differential  algebraic 
equations  (DAEs)  f{x\x,t)  sz  0  where  fx>  is  identically  singular.  The  index  provides 
cme  measure  of  the  singularity  of  a  DAE.  Most  of  the  numerical  analysis  literature 
on  DA^  to  date  has  dealt  with  DAEs  with  indices  no  larger  than  three  because 
of  technical  difficulties  and  because  many  basic  applications  including  constrained 
mechanical  system  have  this  index.  This  paper  will  discuss  several  situations  where 
DAEs  of  index  hi^er  than  three  naturally  occur.  It  will  also  discuss  the  relationship 
between  certain  concepts  in  nonlinear  control  the<»y  such  as  relative  degree  and  zero 
dynamics,  the  index,  and  constrained  mechanical  systems. 

ORAL  PRESENTATIONS  OF  RESEARCH  RESULTS 
1.  Invited  talk  at  Humboldt  University,  Berlin,  June,  1993. 
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2.  at  Conference  on  Scientific  Computation  and  Differential  Equations,  Auckland, 
January,  1993. 

3.  Invited  talk  at  Oberwcdfach  Workslu^  on  Simulation,  June,  1993. 

4.  Invited  talk  at  World  Congreas  of  Nonlinear  Analysts,  Tampa,  1992 

5.  Two  invited  talks  at  1992  SIAM  Summer  Meeting,  Los  Angeles 

6.  Tklk  at  SIAM  Control  Meeting,  Minneapolis,  1992. 

7.  Presented  two  p^>er8  at  the  Symposium  on  Implicit  and  Nonlinear  System  at  the 
Robotics  Institute  of  the  University  of  Texas  at  Arlington,  1992. 

8.  Presented  paper  at  the  31  IEEE  Conference  on  Decision  and  Control,  'Dicson,  1992. 

9.  Gave  talk  at  the  ARO  Workshop  held  in  Raleigh,  1992. 

10.  Invited  talk  at  University  of  Umea,  Sweden,  1991 

11.  Two  invited  talks  at  IMACS  World  Ctmgress  on  Scientific  Computation,  Dublin,  Ire¬ 
land,  1991 

12.  Invited  talk  at  Inst.  Nat.  Res.  Info.  Auto.  (INRIA),  Rocquencourt,  France,  1991 

13.  Invited  talk  at  International  Institute  of  Applied  System  Analysis  (IIASA),  Laxen- 
burg,  Austria,  1991 

14.  Invited  talk  at  Applied  Linear  Algebra  Meeting,  Northern  Illinds  University,  1991 

15.  Eur<q>ean  Control  Conference,  Grenoble,  France,  1991  (panel  discussion) 

16.  Talk  at  29  IEEE  Conference  on  Decision  and  Control,  Hondulu,  1990. 

Stephen  L.  Campbell 
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N(»th  Carolina  State  University 
Raleigh,  NC  27695-8205 
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